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Abstract. We report on the theoretical investigation of Wigner distribution function 
(WDF) reconstruction of the motional quantum state of large molecules in de Broglie 
interference. De Broglie interference of fullerenes and as the like already proves the 
wavelike behaviour of these heavy particles, while we aim to extract more quantitative 
information about the superposition quantum state in motion. We simulate the 
reconstruction of the WDF numerically based on an analytic probability distribution 
and investigate its properties by variation of parameters, which are relevant for the 
experiment. Even though the WDF described in the near-field experiment cannot be 
reconstructed completely, we observe negativity even in the partially reconstructed 
WDF. We further consider incoherent factors to simulate the experimental situation 
such as a finite number of slits, coUimation, and particle-slit van der Waals interaction. 
From this we find experimental conditions to reconstruct the WDF from Talbot 
interference fringes in molecule Talbot-Lau interferometry. 
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1. Introduction 

The reconstruction of a quantum state in general is important to test potential quantum 
systems on their properties such as entanglement, superposition and coherence and 
on their applicability for emerging quantum technologies. Furthermore, massive and 
spatially extended quantum systems such as complex molecules and nanoparticles in 
quantum superposition are the proposed test embodiments for universal boundaries of 
the validity of quantum theory with strong relevance for future nanotechnology [H E] . 
The wave nature of molecular motional states (as massive as ten Cgo molecules) has 
been demonstrated [3] and to make a step further, we now want to characterize the 
quantum state of motion of molecules in Talbot-Lau interferometry. 

The wave function cannot be observed directly, but by using the Wigner distribution 
function (WDF) we have an alternative perspective on quantum dynamics as WDF is 
equivalent to the density matrix |ll|5]. Quantum states have the unique property that 
they can generate negative values of this quasi-probability function. On the first hand 
the negativity of the Wigner function is seen as a proof of the quantum nature of the 
state under consideration. A fully reconstructed Wigner function contains the complete 
information about the measured state. The process to evaluate the state is called phase 
space tomography. 

Phase-space tomography for the Wigner function was pointed out in a general 
context by Bertrand and Bertrand [6] and independently by Vogel and Risken [7], apphed 
in photonics [S] where quantum state tomography is an established experimental tool 
to quantify the quantum state of light pj. It has been recently used to demonstrate 
the state squeezing of atomic Bose condensates PI)]), and to quantify the motional 
quantum state of a trapped Be"*" ion An early experiment to characterize the 

quantum nature of atomic motion in de Broglie interference was the Wigner function 
reconstruction of meta-stable He atoms diffracted at a double slit [12] , as theoretically 
proposed earlier [13]. The Wigner function reconstruction has been also proposed to be 
of use to prove the quantum nature of the superposition of very massive particles [H] and 
even macroscopic opto-mechanical systems [15]. Furthermore, to prove entanglement in 
superconducting qubits, quantum state tomography has been applied [TB]. In more 
technical terms the Wigner function is a quasi-probability distribution of states in 
phase space. In the case of de Broglie interference of molecules the quantum state 
is projected on to the spatial coordinate, which is the spatial number distribution of 
particles after the diffraction grating (see Figjl]). This spatial distribution is needed for 
different rotational angles of the quantum state in phase space. This rotation in phase 
space comes natural from free space propagation of a particle beam after diffraction and 
can be evaluated from measuring the spatial distribution at different distances after the 
grating. The inverse Radon transformation of collected data of density distributions 
will give the Wigner function. In the experiment the spatial coherence is prepared 
by another grating placed in front of the diffraction grating. Then many coherent 
sources constructively contribute to the same interference pattern due to the Lau effect. 
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Figure 1. The setup as used for simulations. The probabihty distribution for 2 
Talbot length z^, the quantum carpet, simulated with finite grating and grating 
opening fraction (slit width / period) of fo = 0.3 is illustrated. The reappearance 
of the grating self-images is the Talbot effect, which occurs if a periodic structure is 
coherently illuminated. 



The multi grating configuration is the so-called Talbot-Lau interferometer (TLI) and is 
explained in detail elsewhere [T71 [18] . 

Here, we discuss the Wigner function reconstruction of matter-waves in the near- 
field Talbot regime with illumination of the grating by a single coherent source. If the 
spatial coherence of the matter wave is high, which means that plane waves are reaching 
the grating, then the here presented Talbot simulations are valid for the Talbot-Lau 
scheme and reconstruction gives the same Wigner function. 

We theoretically perform phase-space tomography of the center of mass motion 
quantum state of massive molecules in a near-field TLI [19j. Near- field means that 
the interference pattern is on the same size scale as the diffraction grating period. In 
contrast, the far-field (Fraunhofer) interference pattern are much larger than the grating 
(for illustration see Hornberger et al. [SO])- Earlier investigations of de Broglie quantum 
states have dealt with far-field pattern reconstruction [12l [13] . From light and matter- 
wave optics it is known that complex diffraction pattern can be expected in this near- 
field regime. Those structures are sometimes called quantum carpets [211 122] • Recently 
optical quantum carpets have been experimentally observed by Case et al. [23]. Talbot 
carpets have been applied for computation as number factorization [211 [2S1 [2S|- 

2. Theoretical Model 

2.1. Phase Space Tomography 

The WDF W{x,p) of a complex signal iplx) is defined by [1]: 

1 fOO 

Wix,p) = ij*{x- x')ij{x + x')e-'^P^ dx', (1) 

TT J-oo 

with momentum p and position x. Throughout the paper we set h=l. Experimentally, 
we measure the spatial intensity pattern which corresponds to the projection of the 
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WDF onto the space coordinate. This is formulated by integration of the WDF over 
the momentum variable p: 



P[x) = / W{x,p)dp. (2) 

J —oo 

If the WDF is rotated by angle 6 in phase space, it becomes 

Wg{x,p) = W{x cos 9 — psm9,x sin 6* + pcos 9), (3) 

and in analogy to Eq. ^ the spatial intensity pattern Pg{x) or the marginal probability 
can be obtained by 

POO 

Pe{x) = / Weix,p)dp. (4) 



This intensity pattern with various angles of rotation can be obtained from diffraction 
pattern like a quantum carpet for a grating with an infinite number of slits. Phase 
space tomography is based on the transformation of Eq. Q resulting in a reconstructed 
WDF: 

1 roo pTT roo 

W(X,P) = — dx' d9Pg(x') rf^|^|eMx.-xcos9-psine)^ /gN 

Atc^ J-oo Jo J-oo 
which is the inverse Radon transformation. Eq. ^ can be written as 

Peix) = X!Pn|^n,e(a;)P, (6) 

n 

with 

^^4x) = ,-1— r dx'e-'^^^'-'^^°''^"^Mx'), (7) 
V Zm sm 9 J-oo 

which represents the fact that the active rotation of the Wigner function can be done by 
fractional Fourier transformation in Fresnel diffraction theory and was derived in |13j . 

2.2. Wigner Function Reconstruction for Free Space Propagation 

A full reconstruction of the WDF requires spatial probability distributions Pe{x) for 
every angle 9 between and tt. The easiest way to rotate the WDF is the free space 
propagation of the particles. The rotation angle then depends on the distance z between 
the diffraction grating and the detector. The free propagation of the diffracted wave 
function is according to Fresnel diffraction theory given by [20] : 

ij{x,z) = r rfxVA("-"')V(a;',0), (8) 

where k = 2tt/\ is equivalent to momentum,with ^ = 1, and A is the de Broglie 
wavelength. We rescale the and get the expression: 

^{±,z) = e-"'^ c/x'e-*'=(^-^V(a;',0). (9) 

S J-oo 

Comparing Eqs. ([T]) and (|9| gives s = ksm(9)/z, for the rescaling of the a;-axis, and 
cot(^) = k/z for the dependency of the rotation angle 9 from the distance z. As a first 
result we find that free space propagation does not lead to a full rotation over vr for finite 
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z. Sufficient rotation to fully reconstruct the WDF can be achieved by a lens. Such 
a lens has been implemented for atomic matter-waves by a Fresnel zone plate [271 l2H] 
or a standing light wave [29]. In [13] it was shown theoretically that by using a lens 
a '7r/2-rotation, and therefore the Fourier transform is accessible for finite values of z. 
Since the realization of such a lens has not been demonstrated for molecular matter- 
waves, we here consider the simplest case of rotation, which is free space propagation. 
However it is possibility to increase the accessible angle of rotation and therefore the 
amount of information about the quantum state by additional symmetry assumptions 



on the investigated WDF [30], which we will use in Sec. (3.1). We will investigate the 
partial reconstruction of WDF in the Talbot regime for B between and tt/2. As a 
figure of merit we use the appearance of negativity of WDF. 

Technically, we numerically reconstructed WDF using the ffitered back-projection 
algorithm [31] for the inverse Radon transformation: 

W{x,p) = dO dx' Pe{x')g{x' — X cos6 — psinO), (10) 

Jo J —Xm 

with 

,1 cos(rcx) rcSin(rcx), , , 

g{x) ^ 2 -- + + ^'' ), 11 

which is approximated from g[x) = dr\r\e'^^^ , where real range of the 

transverse and Tc is a cut-off frequency. 

2. 3. WDF of Ideal Quantum Carpet and Talbot Effect 

The WDF for a quantum superposition state as in the classic double slit arrangement, 
a cat state, is well known. Here, we are interested in the WDF of an ideal quantum 
carpet, the near-field wave diffraction pattern after a grating with many slits. We start 
with a grating with an infinite number of slits, as it gives a nice analytical expression. 
The wave function tplx) after the grating is: 

oo 

tp{x) = tcix)^ 5{x-nd), (12) 

n=— oo 

where d is the grating period, tc{x) is the grating transmission function for a single slit 
{—d/2 < X < d/2), and ® denotes convolution: (/ g){x) = f{x — x')g{x')dx' . An 
infinite train of delta function, a comb function, can be expressed as ^ Z)5^-oo 



With that we rewrite Eq. (12) to be: 



^(x)= E Ane^'-^K (13) 

n=— oo 

where An = d~^ f''^^^i^dxtc{x)e~'^'^'^'^^/'^. The wave function along z is then in form of 
Eq.(g given by: 



ij{x,z)^ J2 A„e-'2."Se-'""'5^. (14) 
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The second exponential term generates a periodic appearance of the identical pattern 
at different positions in z-direction. This is the Talbot effect. Self-images of the grating 
revive at multiple of the Talbot distance zt = 2^. The probability distribution of this 
Talbot effect is shown in FigjlJ We note that self images are also revive at multiple 
of zt/2 with a d/2 displacement in x-direction. It is this d/2 shift which carries non- 
classical information about the motional state and generates negative values of WDF. 
WDF reconstruction will therefore need to include especially such 2;-positions. In the 
following we define the diffracted wave function in units of the Talbot distance, since the 
self-image is repeated at multiple of zt- We note that the number of near-field Talbot 
self-images typically depends on the number of contributing (coherently illuminated) 
grating slits. Experimentally, five Talbot lengths zt are easy to achieve [20] . 

The exact WDF for the infinitely periodic input is obtained from Eq. ([T| as 

oo oo n n' 

W{x,p)= E ^nA:'e^'-(— ')t5(p/27r-^— ). (15) 

n=—oon'=—oo 

The WDF of the (5-comb function, which is the mathematically simplest guess, is 
obtained by substituting m = n + n', and setting A„ = l/c?: 

^ n=-oo m=~oo '^^ 

and plotted in phase space as shown in Fig. [2]^a). The WDF is the sum of 5-function in 
momentum, as seen in Eq. (15). This delta function peaks, however, do not appear in 
experimental physical system where we have a finite number of slits of finite width. A 
grating with finite slit width is formulated with the Fourier coefficient = foSinc{mifo) 
and is plotted in Fig. |2](b). We note the negative peaks (bright spots) in the WDF, 
showing the non-classical character of the state. We observe that the peak width in 
momentum becomes narrower as the number of slits increases. 



3. Simulations: Limits for Finding Negativity in WDF 

Here we simulate the effects of experimental limitations on the quality of the 
reconstructed WDF, which include: 

• Range of rotation angles, 

• Spatial detector resolution, 

• Incoherent source (in terms of collimation of molecule beam, spatial coherence), 

• Finite gratings, 

• Van der Waals interaction (interaction between molecules and material gratings), 

• Visibility. 
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Figure 2. Exact Wigner Distribution Function of (a) (5-function comb wave and (b) 
rectangular wave. Those correspond to the WDF of infinite gratings. The symbols +, 
- in (a) indicate positive, negative delta peaks and the bar in (b) indicates gray scaled 
of WDF in the contour plot. 



We will conclude each section with the feasibility of these limits for the experiment. For 
numerical simulations, the grating period d is set as unit length and 1/d is set as unit 
momentum. All parameter (x, 2;, A, p) become dimensionless by rescaling: x — )■ x/rf, 
z — i- z/d^ A — )■ A/(i, and p — )■ pd. The wavelength is chosen to be A = lO^^c? and 
the open fraction of the gratings is set to /o=0.3. This choice of parameter adapts 
the simulation to molecule interferometry experiments, which are in the centre of our 
interest [IH]. However our results are universal and can easily be tuned to represent 
the same diffraction effects in the Talbot regime for other electromagnetic or matter 



waves. The cut off frequency Tc as defined in Eq. (11) was optimized to = 30 to 
show the structure of WDF without high frequency computational noise. The WDF 
is reconstructed basing on the normalized probability distribution and will be plotted 
as contour maps in all following figures where the gray scale bar indicates the value of 
WDF. The scale bar allow for comparison of the different effects as the intensity scale 
has been normalized for each effect. For some WDF plots we show a cross section of 
the contour plot to visualize the negativity. 

3.1. Rotation Angle 

The WDF is rotated by free space propagation. In the experiment the rotation angle is 
described by the distance z after the grating: 
2n 

z = —tajie. (17) 
A 



Since we use z and A in units of Eq. (17) becomes z/ zt = vrtan^^ which is illustrated 



in Figjsj^a). As a result a rotation of close to 7r/2 can be achieved with several Talbot 
distances zt, although a rotation of exactly 7r/2 corresponds to infinite distance. A 
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Figure 3. (a) Propagation distance z as a function of the rotation angle 9, with 
distance z scaled by Talbot distance zt- (b) Full reconstruction of WDF, and partial 
reconstruction of the WDF with rotation angle between [0,4zt] (c) and [0, zy] (d). 



full reconstruction of the WDF is therefore not feasible, and we now investigate the 
partial reconstruction for limited rotation angles to investigate whether this still leads to 
neg ativity of WDF. In Fig|3| we compare full reconstruction (b), partial reconstruction 
of WDF with rotation angle [0,4^7^] (c), and for [O,^;^] (d). We observe tilting of WDF 
for partial reconstruction, but the position of the peaks and the negativity still remain. 
When partially reconstructing the WDF, negative peaks appear when including the 
displaced self-images, which are truly due to interference. We conclude therefore, these 
negative peaks in the partially reconstructed WDF show non-classical behavior. In the 
following, we chose a rotation angles between and 4zt which corresponds to rotation 
angles between to arctan[4/7r] according to Eq. (17). We note the rotation angle 
depends on the unit length and the corresponding WDF is also expressed as a function 
of the unit length [30] . 
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Figure 4. Partial reconstruction of WDF between to Azt varying resolution of 
X-axis, which corresponds to the spatial resolution of the detector in the experiment, 
(a) Contour plot for dx = O.Old, (b) Cross section plot for WDF at (p/27r)d = O.bd 
for dx = O.Old (red solid line), dx = O.Q5d (green dotted line), and dx = O.ld (blue 
dash-dotted line) where N0 = 100 and /o = 0.3. 



3.2. Resolution in x and z 

In the last section we investigated the range of rotation angles necessary for observing 
negativity in the WDF. Now, our interest is the dependency on the resolution in x- as 
well as in z-direction. The resolution in x-direction corresponds to the spatial resolution 
of the detector. Simulations are shown in Fig. |4| where we compare dx = 0.01c?, 
dx = 0.05d, and dx = O.ld. We conclude that a resolution of ten measurement points 
per grating period d (and possibly even less) are sufficient to reconstruct negativity 
of WDF. This resolution is typically achieved in state-of-the-art molecule interference 
experiments. However structure of the reconstructed WDF becomes a more pronounced 
for higher x-resolution. We chose dx = O.ld for the following simulations. 

We find that an important parameter is the number of rotation angles Ng which 



corresponds to the propagation distance z, as discussed in Sec. |3.1[ On the other hand 
the resolution of the ^-direction itself is not critical and we show simulations in Fig. |5] 
to prove this statement. We keep the the total distance z constant, but change the 
resolution which is the number of scans within z. In Fig. |5](a) the reconstructed WDF 
for Ng=20 within Azt is shown. The structure is very similar compared to Fig. [3|c), 
which is over the same z distance. We then vary the resolution further and show cross 
sections of the contour plot of the WDF at ip/2Tc)d = 0.5 for Ng = 20,50,100 in 
Fig. [sj^b). The periodic negative and positive peaks are almost independent from Ng. 
Therefore, data taken at the self-image planes are sufficient to reconstruct WDF. That 
means, when measuring only the self-image planes we need to measure until 4zt to find 
negativity in the reconstructed WDF. We note, displaced self images taken at integer 
multiple of zt/2 are crucial to generate negative peaks. 
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Figure 5. Partial reconstruction of WDF from to Azt varying the number of rotation 
angles (a) Contour plot for Ng — 20, (b) Cross section plot for WDF at (p/27r)d = 0.5 
for Ng = 100(red solid line), Ng = 50 (green dotted line), and Ng — 20 (blue dash- 
dotted line) where /o = 0.3 and dx = 0.1. 



3.3. Incoherent Source 

Effects of temporal coherence, which correspond to the longitudinal z-velocity selection 
of the molecular beam is taken into account by a reduction of the visibility for a given 



distance z in Sec. (3.6). Here, we discuss the effect of an incoherent source with 
respect to spatial coherence of the molecular matter waves. This corresponds to the 
transverse x-velocity selection by coUimation of the molecular beam. Assuming an 
incident coUimation of angle a (see Fig. ([T])), the wave function is found by averaging 
over all incident angles a: 

I °° 97r (k A-n — r)'^ \ 

ij{x,z) = ^Y. A^exp[z{K + n^x)+z{k-^-^^^^j^)z]'l , (18) 

where k = {kz^ + {ka + n'^x)^)^/'^ and k^ = A;sin(a). The variation of the beam 
coUimation is implemented mathematically by the standard deviation of a Gaussian 
distribution of a = Q.lamax, with amax is the total range of incident coUimation angles. 
The wave function is then averaged between [—amaxl'^ , cxmaxI'A with Gaussian weight. 
The Talbot effect can be observed in the limit a < d/zT, which means that the diffraction 
is dominating coUimation. The probability distribution is plotted in Fig. |6](a) for incident 
angle a^ax = vr x 10~^, which is easy to achieve in the experiment by coUimation at the 
fist grating in a Talbot and Talbot-Lau interferometer. We observe that interference 
contrast is washed out rapidly for increasing z. We reconstruct WDF as shown in 
Fig. |6](b) and it still shows negative peaks between positive peaks. However, we observe 
a significant lost of contrast of WDF for increased coUimation angle a^nax, as shown 
in Fig. ^c) and (d). Since the interference pattern are very sensitive to the spatial 
coherence, the Talbot-Lau interferometer configuration is used in experiments. Here, 
an additional grating placed in front of the diffraction grating increases the spatial 
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Figure 6. Probability distribution (a) and partial reconstruction of the WDF with an 
incoherent source based on data of [0, Azt] with an incident angle of [a]maa; = tt x 10~^ 
(b) and [a]^ax = 2.57r x W-^{c). (d) The WDF at (p/27r)rf = 0.5 for a,nax = (red 
solid line), amax = tt x 10^^ (green dotted line), and 5n x 10^^ (blue dash-dotted line) 
where /o = 0.3, iV^ = 100, da: = O.ld. 

coherence of the matter waves and Talbot images can be well revived even with a 
molecular beam source of relatively poor spatial coherence [19]. 

3.4- Finite Gratings 

The grating was assumed to have an infinite number of slits until now. Experimentally, 
the effective number of involved slits is finite, while sufficiently large to observe the 
Talbot effect. Talbot self-images, however, are washed out after some multiple of 
Talbot distance zt as an effect of the finite number of slits [20j. To study this effect 
quantitatively, we define the wave function for a finite number of slits as: 

n=Ns 

ijj{x) = tc{x) ^ ^ 6{x-nd), (19) 

n=-Ns 
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Figure 7. Probability distribution (a) and partial reconstruction of the WDF (b) for 
finite number of slits, iV,=f for [0, 4zt] where Ne = 100, /o = 0.3, dx = 0.1. 



with 2Ns slits. The wave function after the grating in free space propagation is then 
given by: 



n=Ns 



(20) 



We simulate density distribution and WDF for various number of slits. Contrast in 
the density distribution is still visible in Fig. for as few as ten grating slits Ns 
contributing to diffraction for up to three Talbot distances and we can still reconstruct 
WDF with negativity in Fig. [7]^b). Both, the visibility of positive and negative peaks 
of WDF as well as the distance (number of zt), where near-field pattern can be seen in 
the density distribution, increase as the number of slits increases. 



3.5. Van der Waals Interaction 

Here we investigate the effect of the interaction between the grating and molecules on the 
WDF reconstruction. Dispersion forces, such as van der Waals (vdW) forces in the short 
range limit (on the order of 100 nm between molecule and grating wall) and Casimir- 
Polder (CP) forces in the long range limit (for larger distances than some 100 nm), are 
known to affect the molecule interference pattern [17], if the gratings are realized from 
material structures made of metals (gold, Au) or semiconductors (silicon nitride, SiN^,). 
Dispersion forces will therefore have a large effect on the quantum carpet structure and 
we have to test the respective limits of WDF reconstruction. We implement the grating 
wall molecule interaction as a phase term 0(x) to transmission function t'^{x) = tc{x)(f){x) 



in Eq. (12): 



0(x) = ei^/-°°^'^^^(^'^\ (21) 

with the mass m of the molecule and the VdW potential V = —C^/x^, which we 
exclusively discuss here for simplicity. The 1/x^ scaling of the potential with distance x 
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Figure 8. The partial reconstructed WDF for [0, 4zt] witliout van der Waals 
(vdW) interaction (a) and with vdW interaction (b) where open fraction /o = 0.44, 
dx = 0.1d,Ng = 80. 



is typical for a pointlike particle in front of a surface. The vdW interaction constant C3 
is depending on the dielectric properties of the molecule and the grating over the full 
electromagnetic spectrum. For our simulations we use 6*3=10 meV nm^ for a fuUerene 
Ceo molecule close to a gold surface in agreement with literature JT7] . Fig. [s] shows 
reconstructed WDF (a) without and (b) with vdW interaction and it shows a small 
change in visibility and structure. We note, that density distribution (quantum carpet) 
and WDF with vdW are very similar to simulations with smaller open fraction /q = 0.1 
without vdW. This is in agreement with the known effect that the always attractive 
vdW interaction effectively reduces the slit width of the diffraction grating |18]. The 
reduced visibility in the density distribution and the reduced negativity of WDF in 
Fig. |8] (b) are explained by dephasing. 



3. 6. Visibility 

Another parameter which will be investigated here is the visibility V = {Imax — 
Imin) / {Imax + Imin) , where / is the intensity of the quantum carpet structure illustrated 
by the gray scale in the plots. The visibility is the most sensitive parameter of the pattern 
an is easily influenced by different effects. For instance a low velocity selection of the 
molecular beam resulting in a low longitudinal or temporal coherence of the matter wave 
reduces the visibility of the carpet. The negativity of WDF is therefore also reduced. 



Furthermore as already mentioned in Sec. 3.3, the Talbot effect is experimentally usually 



observed in the Talbot-Lau interferometer. The main difference between Talbot and 
Talbot-Lau from the theoretical point of view is a reduction in interference fringe 
visibility for the TLI. This means our discussion on WDF reconstruction in the Talbot 
regime can be also extended to the Talbot-Lau regime, if we investigate the effect of 
visibility reduction. To include the visibihty in the probability distribution (quantum 



Phase Space Tomography of Matter- Wave Diffraction in the Talbot Regime 14 
(a) (b) 

0.8 1 1 1 1 1 r 

0.6 - 

0.4 - 

^ 

"^-0.2 - 

-0.4 - 

-0.6 - 

-0.8 I ' ' ' ' ^ 

-3-2-1012 

x/d 



Figure 9. (a) Contour plot of partially reconstructed WDF for [0, 4zt] with visibility 
V =0.5, and (b) cross section plot of the WDF at {p/2n)d — 0.5 for V = 1 (red 
solid line), V = 0.75 (green dotted line), and V — 0.5 (blue dash-dotted line) where 
Ng = 100, /o = 0.3, = O.ld. 

carpet) we modify the density distribution P{x.z) to be: Imin + (1 — Imin)P{x, z), with 
Imin = (1 ~ + '^)- We reconstruct WDF and plot in Fig. |9](a), vary the visibihty 

and find that with a visibihty of V=0.5 we can stiU identify negative parts of the WDF 
although the contrast is reduced as shown in Fig. |9](b). The analysis of visibility by 
WDF reconstruction will become a useful tool to investigate decoherence effects and 
mechanisms in molecule quantum optics. 

4. Conclusion 

We have numerically analyzed the partial Wigner distribution function (WDF) 
reconstruction of near-field optical quantum carpets for free space propagation of the 
wave. In our study, we considered all the major experimental inefficiencies. Most 
important, we find that negativity can be observed in the reconstructed Wigner 
distribution function in the Talbot and Talbot-Lau regime under such realistic conditions 
with today's technologies. All investigated parameter and effects keep promise for the 
realization of molecular quantum optics tomography. Important to reconstruct WDF is 
to collect data of quantum carpets at half-integer multiple of Talbot distances. The 
Talbot regime is important for molecule interferometry as the conceptually simpler 
double slit and far-field interference are much harder to perform due to experimental 
difficulties. The here presented tomography of the motional quantum state will 
become an important analytic tool in molecule quantum optics as it gives a handle 
to directly detect the superposition state of centre of mass motion. Furthermore WDF 
reconstruction may be used to investigate wave function dephasing effects such as van 
der Waals or Casimir-Polder interactions of the particles as well as to study decoherence 
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effects which reduce the visibihty V. Further work needs to conduct experimental 
tomography of the discussed quantum state, while the main challenge lasts to generate 
an intense beam of large particles while keeping sufficient coherence, which requires 
high phase-space density. The here described tomographic analysis is also applicable 
to investigate the quantum superposition of the centre of mass motion of even larger 
nanoparticle and cluster as well as for electromagnetic waves such as X-rays [32] . 
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